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Abstract 

We explore the relationship between exclusive and inclusive electromagnetic scattering from 
' ' the pion, focusing on the transition region at intermediate . Combining Drell-Yan data on 

. the leading twist quark distribution in the pion with a model for the resonance region at large 

\ X, we calculate QCD moments of the pion structure function over a range of Q^, and quantify 

^ [ the role of higher twist corrections. Using a parameterization of the pion elastic form factor and 

Qh' phenomenological models for the n ^ p transition form factor, we further test the extent to which 

' local duality may be valid for the pion. 
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I. INTRODUCTION 



The nature of the transition between quark and hadron degrees of freedom in QCD is 
one of the most fundamental problems in strong interaction physics. This transition has 
been extensively explored within nonperturbative models of QCD, which, while retaining 
some of the apposite features of QCD, make simplifying assumptions that allow approximate 
solutions to be found Considerable progress has also been made recently in calculating 
hadronic properties directly from QCD via lattice gauge theory, and much is anticipated from 
this approach in the near future with significant advances in computing power available ■ 
It is clear, however, that while a quantitative description of hadronic structure from first 
principles in QCD is still some time away, phenomenological input will remain crucial in 
guiding our understanding for the foreseeable future. 

Of course, assuming QCD can ultimately describe the physics of hadrons, the transition 
from quarks and gluons to hadrons can be considered trivial in principle from the point 
of view of quark-hadron duality. So long as one has access to a complete set of states, it 
is immaterial whether one calculates physical quantities in terms of elementary quark or 
effective hadron degrees of freedom. In practice, however, truncations are unavoidable, and 
it is precisely the consequences of working with incomplete or truncated basis states that 
allows one to expose the interesting dynamics that drives the quark-hadron transition. 

The duality between quarks and hadrons reveals itself in spectacular fashion through the 
phenomenon of Bloom- Oilman duality in inclusive lepton-nucleon scattering, eA^ —>■ eX. 
Here the inclusive F2 structure function of the nucleon measured in the region dominated by 
low-lying nucleon resonances is observed to follow a global scaling curve describing the high 
energy data, to which the resonance structure function averages 0, 0]- The equivalence of 
the averaged resonance and scaling structure functions in addition appears to hold for each 
prominent resonance region separately, suggesting that the resonance-scaling duality also 
exists to some extent locally. 

The correspondence between exclusive and inclusive observables in electroproduction was 
studied even before the advent of QCD jHSSH, Within QCD, the appearance of 

duality for the moments of structure functions can be related through the operator product 



expansion (OPE) to the size of high twist corrections to the scaling structure function |ll 



which reflect the importance of long-range multi-parton correlations in the hadron 3| . The 
apparent early onset of Bloom-Oilman duality for the proton structure function seen in 
recent Jefferson Lab experiment s |^ indicates the dominance of single-quark scattering to 
rather low momentum transfer [l3|. It is not a priori clear, however, whether this is due 
to an overall suppression of coherent effects in inclusive scattering, or because of fortuitous 
cancellations of possibly large corrections. Indeed, there are some indications from models 
of QCD that the workings of duality may be rather different in the neutron than in the 
proton 3 3|i f*^^ spin-independent and spin-dependent structure functions 



From another direction, one knows from the large N^. limit of QCD that duality is an 
inevitable consequence of quark confinement; in the mesonic sector one can prove (at least in 
1-1-1 dimensions) that an exactly scaling structure function can be constructed from towers of 
infinitesimally narrow mesonic qq resonances |17|]. This proof-of-principle example provides 
a heuristic guide to the appearance of the qualitative features of Bloom-Oilman duality, 
and has been used to motivate more elaborate studies of duality in quark models, even 
though application to the baryon sector is somewhat more involved Oiven that Bloom- 
Oilman duality is empirically established only for baryons (specifically, the proton), while 
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the application of theoretical models is generally more straightforward in the meson sector, a 
natural question to consider is whether, and how, duality manifests itself phenomenologically 
for the simplest qq system in QCD — the pion. 

As the lightest qq bound state, the pion plays a special role in QCD. Indeed, in Nature 
the pion presents itself as somewhat of a dichotomy: on the one hand, its anomalously small 
mass suggests that it should be identified with the pseudo-Goldstone mode of dynamical 
breaking of chiral symmetry in QCD; on the other, it ought to be described equally well 
from the QCD Lagrangian in terms of current quarks, with particularly attractive forces 
acting in the = 0~ channel. The complementarity of these pictures may also reflect, 
loosely-speaking, a kind of duality between the effective, hadronic description based on 
symmetries, and a microscopic description in terms of partons. This duality is effectively 
exploited in calculations of hadron properties via the QCD sum rule method ^3], in which 
results obtained in terms of hadronic variables using dispersion relations are matched with 
those of the OPE using free quarks. 

In this paper we connect a number of these themes in an attempt to further develop and 
elucidate the issue of quark-hadron duality for the pion, focusing in particular on insights 
that can be gained from phenomenological constraints. Specifically, we shall examine the 
possible connections between the structure of the pion as revealed in exclusive scattering, and 
that which is measured in inclusive reactions. The latter can in principle be reconstructed 
given sufficient knowledge of the form factors which parameterize transitions from the ground 
state pion to excited states. We do not attempt this rather challenging task directly; instead, 
we use the tools of the OPE to organize moments of the pion structure function according 
to (matrix elements of) local operators of a given twist. This exercise is possible because 
the structure function of the pion has been determined from Drell-Yan vrA^ scattering data 
at high Q^. Of course, the absence of fixed pion targets means that the structure of the pion 
at low excitation mass W is not known, with the exception of the elastic pion contribution, 
which has been accurately measured for ^ 2 GeV^ in tt"*" electroproduction off the proton. 

To complement the dearth of data on specific it tt* transitions above threshold (but 
below the deep inelastic continuum), we consider a simple model for the pion structure 
function in which the low W spectrum is dominated by the elastic and n —>■ p transitions, 
on top of a continuum which is estimated by evolving the leading twist structure function to 
lower Q^. The discussion at low W is necessarily more qualitative than for the corresponding 
case of the nucleon where ample data exist. However, even within the current limitations, 
this analysis provides an estimate of the possible size of higher twist effects in the pion 
structure function, and the role of the resonance region in deep inelastic scattering (DIS) 
from the pion. Preliminary results for the higher twist corrections have been presented in 
Ref. Here we shall extend that analysis by considering the extent to which local duality 
may be valid in the pion structure function, and possible constraints on the x ^ 1 behavior 
which can be inferred from the elastic channels. 

This study is timely in view of experiments on the pion elastic 23l | and transition 



j2J, |25[ form factors being planned or analyzed at Jefferson Lab, which will probe the 
interplay between soft and hard scattering from the pion and the onset of perturbative 
QCD (pQCD) behavior. Furthermore, recent measurements of the inclusive pion structure 
function via the semi-inclusive charge-exchange reaction, ep —>■ enX, at HERA have yielded 
some unexpected results at low x j26|, and new experiments over a large range of x are being 
planned at Jefferson Lab at lower ^27^. This paper discusses the possible interelations 
between these measurements, in the quest for obtaining a consistent, unified description of 
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the structure of the pion in electromagnetic scattering. 

The structure of this paper is as follows. After briefly reviewing in Section II the defi- 
nitions and kinematics of inclusive lepton scattering from the pion, in Section III we begin 
the discussion by focusing on the special case of elastic scattering. We construct an efficient 
parameterization of the elastic pion form factor in the space-like region consistent with the 

— i> and Q"^ oo constraints. An analysis of moments of the pion structure function 
is presented in Section IV, including the extraction of higher twists and a discussion of the 



role of the resonance region. Some of these results appeared in Ref. [21|. In addition, we 
carefully examine the large x region, which is important for high moments, and compare 
predictions of several models for the leading and higher twist contributions to the pion struc- 
ture function as x — > 1. The relation of the structure function at x ~ 1 with the — *• cxo 
dependence of elastic form factors is discussed in Section V, where we test the hypothesis of 
local Bloom-Gilman duality between the scaling structure function and the exclusive elastic 
and TT ^ p transition contributions. Concluding remarks and a survey of future avenues for 
developments of the issues presented are outlined in Section VI. 



II. DEFINITIONS 

Inclusive scattering of an electron, or any charged lepton, from a pion, evr —>■ eX, is 
described by the pion hadronic tensor, 

W;, = i2nf6\p + q-px)J2{n\J,mX){X\JM\^) , (1) 

X 

where p and q and the pion and virtual photon four-momenta, respectively, and px is the 
momentum of the hadronic final state with invariant mass squared = m\ — q^ + 2m^i/, 
with V the energy transfer in the pion rest frame and the pion mass. The hadronic tensor 
can be parameterized in terms of two structure functions, 

W^, = [-g^, + ^ J ^1 ^ ) + ^^/^ ) ~ -^'i'^ ) —2 ' 

where Wi and W2 are in general functions of two variables, for instance u and g^. In the 
limit as 1/ — 00 and = — — > 00, with x = Q'^/2p ■ q = Q'^/ {W"^ — ml + Q^) fixed, the 
functions Wi and t^H^^ become scale-invariant functions of x, 

m^W^{u,q^) F^{x) , (3a) 
uW^{u,q') ^ F-{x) . (3b) 

Furthermore, in this limit these functions satisfy the Callan-Gross relation, F2 = 2xF^ 
[isj l. Radiative QCD corrections introduce explicit dependence of F^2 the strong cou- 
pling constant, ^^(Q^). While only transversely polarized photons contribute to the F^ 
structure function, F[ oc 0"^, the structure function receives both transverse and longi- 
tudinal contributions, F2 oc ax + ctl, where cxr and ctl are the transverse and longitudinal 
photoabsorption cross sections, respectively. 
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III. PION FORM FACTOR 



The inclusive spectrum begins with the elastic peak at W = m^, or equivalently, x = 
1. Because the pion is spinless, elastic scattering from the pion contributes only to the 
longitudinal cross section, so that the elastic contribution to the structure function 
vanishes. The elastic contribution to the F2 structure function is proportional to the square 
of the elastic pion form factor, 

F^^'''\x = 1,Q2) = 2m^u {F^mY S{W' ~ ml) (4) 

where F^^IQ"^) is the elastic pion form factor. As the most basic observable characterizing the 
composite nature of the lightest bound state in QCD, the elastic form factor of the pion is of 
fundamental importance to our understanding of hadronic structure. In the approximation 
that the pion wave function is dominated by its lowest qq Fock state, the pion form factor 
becomes amenable to rigorous QCD analysis. Indeed, it is well known that the asymptotic 
behavior of the pion form factor is calculable in pQCD 



Fn{Q ) ^ as Q ^ cx) , (5) 

where /tt = 132 MeV is the pion decay constant. Current data on F^r , summarized in Fig. 1, 
indicate that there are large soft contributions still at <2 GeV^ QH^. The low data 
are obtained from scattering pions off atomic electrons 34], while the higher data are 
taken from ^H(e, e'7r+)?T, measurements at CEA/Cornell |j35j, DESY ^] and JLab For 
comparison, the leading order pQCD prediction from Eq. (0) is shown in Fig. 1. Although 
the region of applicability of the pQCD result is a priori unknown, the pion represents the 
best hope of observing the onset of asymptotic behavior experimentally (the corresponding 
pQCD calculation of the nucleon form factors significantly underestimates the data at the 
same Q"^). 

There have been a number of calculations of the elastic pion form factor at low Q^, for 
instance using the QCD sum rule approach [3^. Rather than rely on any specific model, 
however, in this analysis we use empirical data to calculate the elastic contribution to the 
pion structure function. For convenience, and for later use in Sections IV and V, we present a 
simple parameterization of the pion elastic form factor data in the space-like region, which is 
valid over the entire range of currently accessible, and smoothly interpolates between the 
pQCD and photoproduction limits. For the latter, the pion form factor at low can be well 
described in the vector meson dominance hypothesis, in which F-j^iQ"^) ~ 1/(1 + Q'^/m'^p). A 
best fit to the low data using the simple monopole form is shown in Fig. 1 (dashed), with 
a cut-off mass ^0.74 GeV. The monopole fit is not compatible, however, with the behavior 
at high expected from pQCD. Building in the ^ and Q"^ ^ 00 constraints, Eq. 0, 
the available form factor data can be fitted by the form 

F (O^) = ^ f ^ + CiZ + c,Z' \ 

^ 1 + QVml \l + c^Z + c,Z^ + C5Z3 ) ' 

where Z = log(l + Q^/A^), and A is the QCD scale parameter. The form © is similar 
to that proposed in Ref. j39| within a dispersion relation analysis, however, the form there 
uses 2 additional parameters, and takes a rather large value of A ~ 1 GeV. Note that the 
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FIG. 1: Pion form factor as a function of Q^. Shown are the best fits (sohd) using Eq. 0, a 
monopole fit (dashed) with a cut-off mass of 0.74 GeV, and the asymptotic prediction from pQCD 
(dotted). 





Ci C2 C3 C4 


fit I 
fit II 


-0.201 0.020 -0.030 -0.093 
0.100 0.060 0.538 -0.249 



TABLE I: Fit parameters for the pion form factor in Eq. ®, as discussed in the text. 



parameterization © is valid only in the space-like reg ion; for a recent discussion of the 
properties of Ft^{Q'^) in the time-like region see Ref. j40|. 

The best fit parameters C1...4 which give the minimum ^ire given in Table I. The 
parameter C5 is constrained by the pQCD asymptotic limit, C5 = m^(/9o/327r^/^)c2, where 
(3o = 11 — 2Nf/3 (= 9 for the 3-fiavor case) |41[. For the QCD scale parameter we take 
A = 0.25 GeV. For completeness, we offer two parameterizations, which approach the pQCD 
limit (0) differently: in fit I the form factor becomes dominated by hard scattering at 
g2 ^ 100 GeV^, consistent with semi-phenomenological expectations jUj, while in fit II 
around half of the strength of the form factor at this scale still comes from soft contributions. 
These are indicated by the solid lines in Fig. 1. Better quality data are needed to constrain 
F-jriQ'^) at higher (^ 2 GeV^). To this end, there are plans to measure the pion form 
factor at an energy-upgraded Jefferson Lab to ~ 6 GeV^ ^ 
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IV. PION STRUCTURE FUNCTION 



Going from elastic pion scattering {W = m^) to the more general case of inelastic scat- 
tering {W > m^), in this section we analyze the pion structure function, F2, in terms of an 
OPE of its moments in QCD, and obtain an estimate for the size of higher twists corrections 
to the scaling contribution. Following this we discuss the role of higher twists in the pion 
structure function at large x, and compare several models for the x — > 1 behavior of with 
data from Drell-Yan experiments. 



A. Moments 

From the operator product expansion in QCD, moments of the pion structure function, 
defined as 



Jo 



(7) 



can be expanded perturbatively at large as a power series in with coefficients given 

by matrix elements of local operators of a given twist (defined as the mass dimension minus 
the spin of the operator), 

00 / 1 N'^ 

= Y.Al{as{Q'))l^—^j . (8) 

Here the leading twist (twist 2) term corresponds to free quark scattering, and modulo 
perturbative asiQ"^) corrections is responsible for the scaling of the structure functions. The 
higher twist contributions A^^q represent matrix elements of operators involving both quark 
and gluon fields, and are suppressed by additional powers of l/Q^. The higher twist terms 
refiect the strength of nonperturbative QCD effects, such as multi-parton correlations, which 
are associated with confinement. 

Note that the definition of M„(Q^) includes the elastic contribution at x = Q'^/{W^ — 
m\ + Q"^) = 1, where W is the mass of the hadronic final state. Although negligible at high 
Q^, the elastic contribution has been found to be important numerically at intermediate 
for moments of the nucleon structure function ^2^. In the definition ((7j) we use the Cornwall- 
Norton moments rather than the Nachtmann moments, which are expressed in terms of the 
Nachtmann scaling variable, ^ = 2x/ (l + -\/l -|- Ax'^ml/Q'^), that includes effects of the target 
mass. The use of the Cornwall-Norton moments was advocated in Ref. |0| on the grounds 
that it avoids the unphysical region C, > C,{x = 1). Because of the small value of m^, the 
difference between the variables x and ^, and therefore between the x- and ^-moments, is 
negligible for the pion. 

The seminal analysis of De Riijula et al. ^ (see also Ref. ji^]) demonstrated that 



the onset of quark-hadron duality is governed directly by the size of the higher twist matrix 
elements. In particular, duality implies the existence of a region in the (n, Q^) space in which 
the moments of the structure function are dominated by low mass resonances, and where the 
higher twist contributions are neither dominant nor negligible. For the case of the proton 
F2 structure function, even though there are large contributions from the resonance region, 
conventionally defined as < 2 GeV, to the n = 2 moment (~ 70% at = 1 GeV^), the 
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higher twists contribute only around 10-20% to the cross section at the same [20|- The 
question we wish to address here is whether there exists an analogous region for the pion, 
where the resonance contributions are important, but higher twist effects are small enough 
for duality to be observed. 

Of course, the distinction between the resonance region and the deep inelastic continuum 
is in practice somewhat arbitrary. In the large A^^^ limit of QCD, for instance, the final 
state in DIS from the pion is populated by infinitely narrow resonances even in the Bjorken 
limit, while the structure function calculated at the quark level produces a smooth, scaling 
function Empirically, the spectrum of the excited states of the pion is expected to 



be rather smooth sufficiently above the p mass, for W ^ 1 GeV. Resonance excitation of 
heavier mesons are not expected to be easily discernible from the DIS continuum — the oi 
meson, for instance, at a mass W ~ 1.3 GeV, has a rather broad width (~ 350 — 500 MeV) 



Moments of the pion structure function can also be calculated directly via lattice QCD, 
and first simulations of the leading as well as some specific higher twist contributions have 
been performed jHf. Although the detailed x dependence, especially at large x (see next 
section) requires knowledge of high moments 3| , considerable information on the shape of 
the valence distribution can already be extracted from just the lowest 3 or 4 moments jiij ]. 
Calculations of a further 2 or 3 moments may be sufficient to allow both the valence and 
sea distributions to be extracted function of x. 

Measurements of the pion structure function have been made using the Drell-Yan process 
|46l 47, lil, 0, S^l in t^N scattering, covering a large range of x, 0.2 < x < 1, and for 



typically > 20 GeV^. It has also been extracted from the semi-inclusive DIS data at 
HERA for very low x and high W j2^ (see also Ref. jlj). However, there are no data on 
F2 at low W , in the region where mesonic resonances would dominate the cross section. 
The spectrum could in principle be reconstructed by observing low t neutrons produced in 
the semi-inclusive charge-exchange reaction, ep — > enX, where t is the momentum transfer 
squared between the proton and neutron, and extrapolating to the pion pole to ensure vr 
exchange dominance. In the meantime, to obtain a quantitative estimate of the importance 
of the resonance region, we model the pion spectrum at low W in terms of the elastic 
and p pole contributions, on top of the DIS continuum evolved down from the higher Q"^ 
region, as outlined in Ref. I21I . The leading twist structure function can be reconstructed 
from parameterizations (s^. l53l Is^ of quark distributions in the pion obtained from global 
analyses of the pion Drell-Yan data. Unless otherwise stated, in this work we use the low 
fit from Ref. [H^l, which gives the leading twist parton distributions in the pion for 
> 0.25 GeV^ (our conclusions are not sensitive to the use of other parameterizations 
E3|). For the elastic contribution we use the parameterization in Eq. © [fit I]. 
The contribution of the p meson is described by the tt ^ p transition form factor, 
F'npiQ'^), which is normalized such that -F^p(O) = 1, and is expected to fall as l/Q^ at large 
(compared with l/Q^ for F^((5^)). Since there is no empirical information on F^p((5^), 
we consider several models in the literature, based on a relativistic Bethe-Salpeter vertex 
function jH^, a covariant Dyson- Schwinger approach [H^, and light-cone QCD sum rules 




57| . These represent a sizable range (~ 100%) in the magnitude of F^p(Q ) over the region 



of covered in this analysis, with the calculation of Ref. [5|^ giving a somewhat smaller 
result than those in Refs. [ssl . Ejj l. The spread in these predictions can be viewed as an 
indicator of the uncertainty in this contribution. The vr — > p transition form factor can be 
extracted, for instance, from p electroproduction data off the proton, ep — > epp^ 0, |25| . 
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FIG. 2: Contributions to moments of the pion structure function from the resonance region, W < 
Wrcs = 1 GeV, relative to the total |21 1 . 



by reconstructing the decay of the p into two pions. It also forms an important input into 
the calculation of meson-exchange current contributions to deuteron form factors at large 

The contributions from the "resonance region", W < Pl^res = 1 GeV, to the moments of 
the pion structure function, 

MriQ') = f dxx^-'F^{x,Q') , (9) 

J Xi-cs 

are plotted in Fig. 2 as a ratio to the total moment, for n = 2, ■ • ■ , 10. The integration in 
Mn^iQ"^) is from a;i.es = Q'^/iW^cs -fnl + Q"^) to the elastic point, a; = 1. The low W region 
contributes as much as 50% at = 2 GeV^ to the total n = 2 moment, decreasing to 
^ 1% for ^ 10 GeV^ ji^]. Higher moments are more sensitive to the large x region, and 
subsequently receive larger contributions from low W. The n = 10 moment, for example, 
is almost completely saturated by the resonance region at Q"^ = 2 GeV^, and even at 

= 10 GeV^ still receives some 40% of its strength from W < 1 GeV even at = 10 GeV^. 

The relatively large magnitude of the resonance contributions suggests that higher twist 
effects play a more important role in the moments of the pion structure function than for 
the case of the nucleon. In Fig. 3 the lowest {n = 2) moment of is displayed, together 
with its various contributions. The leading twist component, 

M^^m = f dxx--^Fl^^{x,Q') , (10) 

^0 
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FIG. 3: Lowest (n = 2) moment of the pion structure function. The leading twist (sohd) and 
elastic (dashed) contributions are shown, and the shaded region represents the total moment using 
different models for the tt — > p transition 21 1 . 



is expressed (at leading order in ^^(Q^)) in terms of the twist-2 quark distributions in the 
pion, 

Flisri^.Q^) = J]eJxg-(x,Q2) , (11) 
g 

where the valence part of q'^ is normalized such that / dx q^g^i{x, Q^) = 1. The leading twist 
contribution is dominant at > 5 GeV^, while the deviation of the total moment from the 
leading twist at lower indicates the increasingly important role played by higher twists 
there. While negligible beyond ^ 4 GeV^, the elastic contribution is as large as the 
leading twist already at ~ 1 GeV^. The tt —>■ p contribution is more uncertain, and the 
band in Fig. 3 represents the total moment calculated using different models jH, E^l 
of Ft^pIQ"^). However, while the current uncertainty in this contribution is conservatively 
taken to be ~ 100%, doubling this would lead to a modest increase of the band in Fig. 3. 
Uncertainty from poor knowledge of the leading twist distributions at small x ,5^ is 
not expected to be large. 

To extract the higher twist part of the moments, one needs to subtract the leading twist 
contribution in Eq. (jlO|) from the total moments M„(Q^), 

M^^m = M^m - M^^m - u^^m , (12) 

where M™{Q'^) arises from target mass corrections. Because the target mass correction, 
which is formally of leading twist, is proportional to m^/Q^, its contribution will only be 
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(GeV^) 



FIG. 4: Higher twist contribution to the n = 2 moment of the pion structure function, as a ratio to 
the total moment. The band indicates the uncertainty due to the model dependence of the vr — > p 
transition form factor 12 ill. 



felt when 



which is far from the region where the twist expansion is expected 



to be vahd. In principle nonperturbative effects can mix higher twist with higher order 
effects in a^-, rendering the formal separation of the two problematic f£L, SlL, Indeed, 
the perturbative expansion itself may not even be convergent. However, by restricting the 
kinematics to the region of in which the 1/Q^ term is significantly larger than the next 
order correction in a^, the ambiguity in defining the higher twist terms can be neglected 
(2ol |. In Fig. 4 the higher twist contribution to the n = 2 moment is displayed as a function 
of Q^. The band again represents an estimate of the uncertainty in the n p transition 
form factor, as in Fig. 3. At = 1 GeV^ the higher twist contribution is as large as the 
leading twist, decreasing to ~ 1/3 at = 2 GeV^, and vanishes rapidly for ^ 5 GeV^. 



As observed in Ref. 



21 



1 GeV^ appears 
structure functions. 



the size of the higher twist contribution at ( 
larger than that found in similar analyses of the proton F2 |20.] and gi \Qc 
This can be qualitatively understood in terms of the intrinsic transverse momentum of 
quarks in the hadron, {k'^), which typically sets the scale of the higher twist effects. Since 
the transverse momentum is roughly given by the inverse size of the hadron, (k^) ~ l/R"^, 
the smaller confinement radius of the pion means that the average (k^) of quarks in the 
pion will be larger than that in the nucleon. Therefore the magnitude of higher twists in 
is expected to be somewhat larger (0(50%)) than in . The E615 Collaboration indeed 
finds the value {k'^) = 0.8 ± 0.3 GeV^, within the higher twist model of Ref. (g^I- The 
experimental value is obtained by analyzing the a: — >■ 1 dependence of the measured p^ ii~ 
pairs produced in vriV collisions, and the angular distribution at large x. We discuss this in 
more detail below. 
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B. X — > 1 Behavior 



The X ^ 1 behavior of structure functions is important for several reasons. As discussed 
in the previous section, higher moments of receive increasingly large contributions from 
the large-x region, so that a reliable extraction of higher twists from data requires an accurate 
determination of quark distributions at x ~ 1. In addition, since the x ~ 1 region is 
dominated by the lowest qq Fock state component of the pion light-cone wave function, in 
which the interacting quark carries most of the momentum of the pion, the behavior of the 
structure function at x — ^ 1 is expected to be correlated with that of the elastic form factor 
at — > oo. In this section we review various predictions for in the limit as x — 1, and 
relate these to the effects of higher twists discussed in the previous section on the x — > 1 
behavior of the structure function. 

Working within a field theoretic parton model framework which predates QCD, Drell and 
Yan ^ and West showed that if the asymptotic behavior of the form factor is (l/Q^)", 
then the structure function should behave as (1 — x)^"~^ as x — 1. This is referred to as 
the Drell- Yan- West (DYW) relation. Simple application to the case of the pion, in which 
the elastic form factor behaves as 1/Q^ at large Q^, leads to the prediction 

F^ix^l) ~ (1-x) . (13) 



This behavior is also predicted in the model of Ref. 

A dynamical basis for the exclusive-inclusive relation was provided with the advent of 
QCD. By observing that the interacting quark at large x is far off its mass shell, Farrar and 
Jackson JgBI derived the x — *■ 1 behavior of the structure function at x — > 1 by considering 
perturbative one gluon exchange between the q and q constituents in the lowest Fock state 
component of the pion wave function. They found a characteristic ~ (1 — x)^ dependence 
for the transverse part of F2, in apparent contradiction with the naive DYW relation (the 
breakdown of the DYW relation for spinless hadrons was discussed earlier by Landshoff and 
Polkinghorne j^). The longitudinal cross section was found to scale like l/Q"^ relative to the 
transverse jo^. Using so-called 'softened' field theory 6^, in which the pion-quark vertex 



function is described by a Bethe-Salpeter type equation, Ezawa [l3| found a similar (1 — x)^ 
behavior. 

Gunion et al. (g^I later generalized the gluon exchange description to include subleading 
l/Q^ corrections for both the F2 and longitudinal Fl structure functions at x — >^ 1, 

F-(x) ~ S,{l-xf + ^, (14) 

m^) - ^ , (15) 

where the constants S2, T2 and 5*^ are determined phenomenologically. More generally, 
according to the pQCD 'counting rules' js^, the leading components for any hadron with 
n spectator (non-interacting) partons were found to behave as (1 — 3;)2n-i+2|A5z| -^^ 
the X —>■ 1 limit, where ASz is the difference between the helicities of the hadron and the 
interacting quark. More recently, other nonperturbative models have been used to calculate 
the pion structure function |68||, however, because of difficulties associated with incorporating 
high momentum components of the wave function, these may not be reliable in the x ~ 1 
region. 
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FIG. 5: Valence quark distribution in the pion extracted from the FNAL E615 Dreh-Yan experiment 
[5fll |. fitted with leading twist (dashed) and leading + higher twist (solid) contributions, as in 
Eq. The functional forms (1 — x) and (1 — x)^ (dotted) are shown for comparison. 



The predicted x — > 1 behavior of the pion structure function can be tested by comparing 
with Drell-Yan data. The x dependence of the pion quark distributions has been measured 
in Drell-Yan ix^ ijT pair production in ttN collisions (in practice, it A) at BNL [4^], CERN 

50l |. The data for q^i^x) = u'^^ (x) = d'^^ (x) from the most 
are shown in Fig. 5 for 4.05 < m^^ < 8.55 GeV, where m^^ 
pair. The scale dependence within this region was found 



M \M and at Fermilab 



recent Fermilab experiment [5C 
is the invariant mass of the 
to be small [50]. The data were fitted using the form 



q-{x,Q' 



Nx^l-xf + ^f 



2 x^ 
9Q^ 



(16) 



where is a constant, fixed by normalization, and the scale is identified with the 
the dimuon mass squared, m^^. The form parameterizes both leading and higher twist 
effects. Including corrections from evolution, the best fit value for the exponent governing 
the X ^ 1 behavior was found to be 5 ~ 1.21-1.30 [H^l, consistent with the findings of the 
earlier CERN experiments 4^,5 49j- The result of the leading twist E615 fit with b = 1.27 is 
shown in Fig. 5 (dashed). The forms (l — x) 

MM and {1-xf [iiHIellellei are also 

shown for comparison at large x (dotted). The data clearly favor a shape closer to (1 — x), 
rather than the (1 — x)"^ shape implied by the counting rules [30']. 

It has been suggested [6^ that higher twist effects in the pion structure function could 
obscure the true leading twist behavior. The higher twist coefficient T2 was calculated in a 
pQCD-inspired model by Berger and Brodsky [64] in terms of the intrinsic quark momentum 
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in the pion, 



To 



(17) 



that the higher twist contribution may 
X ] as oo, and mimic the observed 



Since it is independent of x, it can be argued [& 
in fact dominate the scahng term at fixed Q^{1 
(1 — a;) dependence if {k^) ~ 1 GeV^. Conway et al. [5i| subsequently performed an analysis 
of the E615 data by fitting also the term 7 in Eq. (fTI)|) . The extracted value of b was found 
to be largely independent of the value of 7 chosen. To investigate whether the quadratic 
term may be masked by an additional component not included in the model [H^l, Conway 
et al. searched for a nonzero intercept of at x = 1. The fit with 7 = 0.83 GeV^ was 
found to be only marginally better than that with 7 = (the significance being 2.5 standard 
deviations), although the fit at x ~ 1 was also sensitive to the input nucleon sea distributions 
in the analysis of the Drell-Yan data. The leading + higher twist fit with 7 = 0.83 GeV^ is 
shown in Fig. 5 (solid curve). The effect on the overall fit is indeed quite marginal, although 
at very large x (^ 0.9) the differences between this and the pure leading twist fit are more 
apparent. 

Mueller j69| has pointed out that Sudakov effects, which introduce terms like 
as(Q^) ln^(l/(l — x)) into the x — > 1 analysis, may invalidate the usual renormalization 
group analysis of DIS at large x. Including power and double logarithmic corrections, one 
finds that the x — > 1 behavior of in this case becomes 



1^1 — x) exp 



In- 



1 — X 



InlnQ^ 



ln^(l/(l-x)) 
21nQ2 



-In- 



1 — X 



In In 



1 



1 — X 



(18) 



When ln(l/(l — x)) = 0{\nQ^ / InlnQ'^) higher twist terms compete with the leading twist, 
and the dominant contribution is then from the longitudinal structure function, which be- 
haves as Fl ~ (1/Q2) ln((5^(l — x)). Taking this criterion literally, for ~ 1 GeV^ this 
would occur at x ~ 0.93, while for Q'^ ~ 100 GeV^ the higher twists would be expected 
to dominate at x ^ 0.97. A cautionary note regarding Eq. (fTH|) . however, is that single 
logarithmic effects have not been included in the analysis, and their effects on Eq. ()18j] are 
unclear. Further discussion of these effects can be found in Refs. jH^ |6^. Carlson and 
Mukhopadhyay [t^ have also studied the effects of radiative corrections on the x — > 1 be- 
havior of the structure function, and the appearance of higher twists in the low-iy region. 
In particular, the scale dependence of the (1 — x) exponent was found to be (1 — x)'^^'^^'^^^'^ , 
with c calculable perturbatively. The dependence of the pion structure function at x ~ 1 
clearly deserves further study. 

A cleaner signature of high twist effects at large x comes from the angular distribution 
of dimuon pairs produced in Drell-Yan collisions. The angular dependence of the Drell-Yan 
cross section is given by [fH (see also Ref. Q|) 



oc 1 + A cos^ 6* + /i sin 29 cos H — sin^ 6 cos 
ail 2 



(19) 



where the angles 6 and are defined in the U~^U~ rest system, and A, fi and u are functions of 
the kinematic variables. In the model of Ref. |6^, the leading twist (1— x)^ term is associated 
with a (1-l-cos^ 9) dependence, while the higher twist {k^)/Q^ term has a characteristic sin^ 9 
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dependence. In particular, the transverse cross section corresponds to A = 1, while deviation 
from a pure (1 + cos^ 6) dependence would indicate the presence of longitudinal or higher 
twist contributions. The data [sOl are consistent with A = 1 for x ^ 0.6, while the larger-x 
data show clear deviations from pure transverse scattering, suggesting the presence of higher 
twist contributions at these x values. With the fitted value of P (~ 1.2-1.3), the measured 
X dependence of A could be accommodated with (k^) ^ 0.8 GeV^. Using the value (3 = 2 
predicted by the pQCD counting rules, the observed A values could be made to fit the data 
by requiring that (fcfn) ~ 0.1 GeV^. However, in addition to being much smaller than the 
value (/cf.) ~ 1 GeV^ suggested in Ref. this scenario is disfavored by a direct comparison 
with the X dependence of g'^(x), as discussed above. 

While the values of the quark intrinsic transverse momentum extracted from the Drell- 
Yan data are consistent with the size of higher twist effects observed in Sec. Ill, there does 
appear to be a clear conflict between the counting rule predictions for the x — * 1 behavior 
of F2 and the empirical x dependence. Several reasons could account for this discrepancy. 
Even higher twist effects, beyond those of twist-4 parameterized in Eq. (jl6|) . could be present 
and obfuscate an underlying (1 — x)^ leading twist behavior. This appears unlikely, however, 
given the relatively large values (Q^ ~ 20 GeV^) at which the data are sampled, and the 
rapid fall off of the higher twist contributions to the moments observed in Sec. IV. A. 

On the other hand, as alluded to above, the extraction of the pion structure function 
requires as input the parton distributions in the nucleon. Since the bulk of the data for 
X > 0.5 corresponds to a nucleon light-cone momentum fraction x^ ~ 0.05-0.1, errors may 
be introduced into the analysis through poor knowledge of the sea quark, or (at higher 
order) gluon, distributions in the nucleon. Furthermore, because the data are taken on 
nuclear targets (e.g. tungsten for the E615 experiment), nuclear effects may give rise to 
corrections to the nucleon quark distributions, especially in the region xa? ~ 0.05, where 
nuclear shadowing is known to play an important role [t^. The effects of using more 
modern nucleon parton distributions, and including nuclear corrections in the analysis, are 
currently being investigated [t^ . 

It may also be that the asymptotic behavior does not set in until x is very close to 1, and 
that the functional form (jl6j) is simply too restrictive to adequately reflect this behavior, 
in which case a more sophisticated parameterization would be required. Further, nonobser- 
vation of the predicted counting rule behavior may not necessarily imply a breakdown of 
pQCD. The derivation of the counting rules for large-x structure functions from Feynman 
diagrams in terms of hard gluon exchanges between quarks involves an infrared cut-off mass 
parameter, m, which regulates the integrals when kr ^ (gtI . ItsI I . Although an analysis 
based on pQCD should be valid also for m = 0, the counting rule results are sensitive to the 
parameter, m, and comparison with phenomenology requires a nonzero value j67l |. 

Regardless of the ultimate x ^ 1 behavior of extracted from data, it is instructive 
to examine whether the asymptotic inclusive-exclusive relations between the pion structure 
function and the pion elastic and transition form factors at large can provide additional 
constraints. In the next section we use local quark-hadron duality to study these relations 
in more detail. 



V. LOCAL QUARK-HADRON DUALITY 

There has been a revival of interest recently in the phenomenon of Bloom-Gilman duality 
in electron-nucleon scattering. This has been stimulated largely by recent high precision 
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measurements |^ at Jefferson Lab of the F2 structure function of the proton, which demon- 
strated that duahty works remarkably well for each of the prominent low-lying resonance 



regions, including the elastic |76|, |77|, as well as for the integrated structure function, to 
rather low values of Q^. Ongoing and planned future studies will focus on duality in other 
structure functions, such as gi J8] and [t^, and for hadrons other than the proton. 

While the existence of loc al q uark-hadron duality appears inevitable in QCD at asymp- 
totically large momenta 0, Isoj . it is not a priori clear that it should work at finite Q^. 
Indeed, there are reasons why at low it should not work at all and its appearance 



in principle be due to accidental cancellations (due to quark charges in the proton 



8l|, for instance) of possibly large higher twist effects. A systematic study of local 
duality for other hadrons, such as the pion, is therefore crucial to revealing the true origin 
of this phenomenon. 

Shortly after the original observations of Bloom-Gilman duality for the proton gen- 
eralizations to the case of the pion were explored. By extending the finite-energy sum rules 
js^l devised for the proton duality studies, Moffat and Snell derived a local duality sum rule 
relating the elastic pion form factor with the scaling structure function of the pion j83l |. 

[FAQ^)f ~ / duo vW^{uj) , (20) 



^max 



where = F2 is a function of the scaling variable uj = 1/x. The upper limit cUi, 

1 + (VF'max ~ ™l)/Q'^ was set in Ref. Q by W^max ~ 1-3 GeV, in order for the integration 
region to include most of the effect of the hadron pole, and not too much contribution 
from higher resonances 18311 . To test the validity of the finite-energy sum rule relation 
(1201), Moffat and Snell j83l], and later Mahapatra [s^l, constructed Regge-based models of 
the pion structure functi on ( their analyses predated the Drell-Yan pion structure function 
measurements jiil 47 . isl . 49l 5^) to compare with the then available pion form factor data. 



The existence of Drell-Yan data on F2 now allows one to test this relation quantitatively 
using on/y phenomenological input. Using parameterizations of the -F;J^(x) data from Ref. [50| 
(see Fig. 5), the resulting form factor -F^(Q^) extracted from Eq. ()2U|) is shown in Fig. 6. The 
agreement appears remarkably good. On the other hand, the magnitude of the form factor 
depends somewhat on the precise value chosen for W^max, so the agreement in Fig. 6 should 
not be taken too literally. Nevertheless, the shape of the form factor is determined by the x 
dependence of the structure function at large x. In particular, while a (1 — x) behavior leads 
to a similar dependence to that for the E615 fit, assuming a (1 — x)^ behavior gives a 
form factor which drops more rapidly with Q^. This simply reflects the kinematic constraint 
(1 - l/u) ~ l/g2 at fixed W. 

Although the apparent phenomenological success of the local duality relation (plj) is 
alluring, there are theoretical reasons why its foundations may be questioned. In fact, the 
workings of local duality for the pion are even more intriguing than for the nucleon. Because 
it has spin 0, elastic scattering from the pion contributes only to the longitudinal cross section 
{F^{x = 1, Q2) = 0). On the other hand, the spin 1/2 nature of quarks guarantees that the 
deep inelastic structure function of the pion is dominated at large by the transverse cross 



section 



111 |65| . Taken at face value, the relation ^J]\i would suggest a nontrivial duality 
relation between longitudinal and transverse cross sections. Whether local duality holds 
individually for longitudinal and transverse cross sections, or for their sum, is currently 
being investigated experimentally. Indications from proton data are that indeed some sort 
of duality holds for both the transverse and longitudinal structure functions of the proton 
individually Hlli. 
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FIG. 6: Local duality prediction for the pion form factor, using phenomenological pion structure 
function input from the FNAL E615 Drell-Yan experiment ,501] (solid), and the forms F2{x) ~ 
(1 — x) and (1 — x)^ (dashed). The asymptotic leading order pQCD prediction (dotted) is shown 
for reference. 



While the elastic form factor of the pion is purely longitudinal, the n ^ p transition is 
purely transverse. It has been suggested jlH that the average of the pion elastic and n ^ p 
transition form factors may instead dual the deep inelastic pion structure function at x ~ 1. 
If we take the simple model used in Sec. IV for the low W part of the pion structure function, 
in which the inclusive pion spectrum at ^ 1 GeV is dominated by the vr — vr and vr — p 
transitions, we can estimate the degree to which such a duality may be valid. Generalizing 
Eq. fj20|) to include the lowest-lying longitudinal and transverse contributions to the structure 
function, one can replace the left hand side of with [Ft,{Q^)]'^ + ct;p[F^p((5^)]^, where 
ujp = l + {ml-ml)/Q\ 

The sum of the lowest two 'resonance' contributions (elastic + p) to the generalized finite- 
energy sum rule is shown in Fig. 7 as a ratio to the corresponding leading twist DIS structure 
function over a similar range of W. The upper and lower sets of curves envelop different 
models 0,0 113 of F,p(Q2), which can be seen as an indicator of the current uncertainty 
in the calculation. Integrating to W^max = 1 GeV, the resonance/DIS ratio at ~ 2 GeV^ 
is ~ 50 ± 30% above unity, and is consistent with unity for ~ 4-6 GeV^ (solid curves). 
As a test of the sensitivity of the results to the value of VFmax, the resonance/DIS ratio is 
also shown for VFmax = 1-3 GeV (dotted curves). In this case the agreement is better for 

~ 1-3 GeV^ with the ratio being ~ 30 ± 20% below unity for ~ 4-6 GeV^. 

Given the simple nature of the model used for the excitation spectrum, and the poor 
knowledge of the n p transition form factor, as well as of the pion elastic form factor 
beyond 2 GeV^, the comparison can only be viewed as qualitative. However, the 
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(GeV^) 

FIG. 7: Ratio of the pion resonance (elastic + vr — > /? transition) contributions relative to the 
DIS continuum, for different values of VFmax- The two sets of upper and lower curves reflect the 
uncertainties in the tt — > p transition form factor. 

agreement between the DIS and resonance contributions appears promising. Clearly, data 
on the inclusive tt spectrum at low W would be invaluable for testing the local duality 
hypothesis more quantitatively. In addition, measurement of the individual transverse and 
longitudinal cross sections of the pion, using Rosenbluth separation techniques, would allow 
duality to be tested separately for the longitudinal and transverse structure functions of the 
pion. 

VI. CONCLUSION 

Understanding the structure of the pion represents a fundamental challenge in QCD. High 
energy scattering experiments reveal its quark and gluon substructure, while at low energies 

its role as a Goldstone boson mode associated with chiral symmetry breaking in QCD is 
essential in describing the long-range structure and interactions of hadrons. We have sought 
to elucidate the structure of the pion by considering its response to electromagnetic probes, 
focusing in particular on the connection between inclusive and exclusive channels. 

The relation between the pion structure function and the pion elastic and transition form 
factors has been studied in the context of quark-hadron duality. Moments of the pion struc- 
ture function have been evaluated, and the role of the resonance region studied, assuming 
that the low W resonant spectrum is dominated by the elastic and tt ^ p transitions. The 
contribution of the resonance region iyV < 1 GeV) to the lowest moment of is ~ 50% 
at «i 2 GeV^, and only falls below 10% for ^ 5 GeV^. The elastic component, while 
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negligible for ~ 3 GeV^, is comparable to the leading twist contribution at ~ 1 GeV^. 
Combined, this means that the higher twist corrections to the n = 2 moment are ~ 50% at 

= 1 GeV^, ~ 30% at = 2 GeV^, and only become insignificant beyond Q"^ ^ 6 GeV^. 

Uncertainties on these estimates are mainly due to the poor knowledge of the inclusive 
pion spectrum at low W, which limits the extent to which duality in the pion can be tested 
quantitatively. Only the elastic form factor has been accurately measured to ~ 2 GeV^, 
although at larger it is poorly constrained. The inclusive pion spectrum can be extracted 
from data from the semi-inclusive charge-exchange reaction, ep — >■ enX, at low t, for instance 
at Jefferson Lab 25j . This could also allow one to determine the individual exclusive channels 
at low W. In addition, a Rosenbluth separation would allow the transverse and longitudinal 
structure functions to be extracted. 

Within the current uncertainties, the higher twist effects in the pion appear lar ger than 
the analogous corrections extracted from moments of the nucleon structure functions^2ol. l63l| . 
This can be generically understood in terms of the larger intrinsic transverse momentum 
of quarks, which governs the scale of the l/Q"^ corrections, in the pion than in nucleon, 
associated with the smaller pion confinement radius. The implication is that duality would 
therefore be expected to set in later (at larger Q^) for the pion than for the nucleon. 

Higher twist effects have also been observed in the pion structure function at large x 
by the E615 Collaboration at Fermilab [s^. The x dependence and angular distribution of 
pairs produced in vrA^ collisions at a: ~ 1 suggests a value {k"^) = 0.8 ± 0.3 GeV^, 
which is larger than the typical quark transverse momentum in the nucleon ((9(500 MeV)). 
On the other hand, the measured x dependence appears to be harder than that predicted 
by counting rules jSOj or models based on perturbative one gluon exchange [11 0, 0, > 
favoring a (1 — x) shape over a (1 — x)^ dependence. A reanalysis j?^ of the Drell-Yan 
data to take into account nuclear corrections and updated sea quark distributions in the 
nucleon, which are used as input into the analysis, is necessary for a definitive assessment of 
the validity of the various approaches. Additional modification of the x ^ 1 behavior due 
to Sudakov-like effects may also need to be considered before drawing final conclusions 
about the implications of the observed a; — *• 1 dependence. There are also plans to measure 
F2 in semi-inclusive reactions over a range of x at Jefferson Lab |i27,J to confirm the Drell-Yan 
and semi-inclusive HERA measurements, which should allow a more thorough exploration 
of the higher twist effects at lower Q^. 

The specific x —>■ 1 behavior of the pion structure function has consequences for the 
dependence of the elastic pion form factor, if one assumes the validity of local quark- 
hadron duality for the pion. In particular, using parameterizations of the Drell-Yan structure 
function data, the existing data on F,r(Q^) can be fitted if the upper limit of the integration 
region above the elastic peak extends to Wmax ~ 1-3 GeV. Analogous fits with a (1 — x)^ 
shape fall off too rapidly with and do not fit Fj^{Q'^) as well. 

On the other hand, there may be limitations of the extent to which local duality can 
hold for the pion, as such duality implies a nontrivial relationship between the longitudinal 
and transverse cross sections. It may in fact be more appropriate to examine whether the 
sum of the longitudinal (elastic) and transverse (vr — p transition) contributions duals the 
DIS structure function at low W . Using phenomenological models for the ir ^ p form 
factor, our estimate for the sum of the lowest-lying resonant contributions is in qualitative 
agreement with the corresponding scaling contribution in the same W interval. However, 
empirical information on the strength and dependence of F-i^p^Q"^) is necessary for a more 
quantitative test. The vr — > p transition form factor can in practice be extracted from p 
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electroproduction data [24]. At larger Q^, the it p transition is expected to be suppressed 
relative to the elastic contribution, and to test the local duality here will require a more 
accurate determination of F-n-fQ^). The pion form factor F^^IQ"^) will soon be measured to 

= 2.5 GeV^ at Jefferson Lab [22| in vr"*" electroproduction from the proton, and possibly 
to = 6 GeV^ with an energy upgraded facility |23l |. 

Finally, this analysis can be easily extended to the strangeness sector, to study the duality 
between the form factor and structure function of the kaon. Data from the Drell-Yan reaction 
in /^"-nucleus collisions [s^ indicate that the quark distribution in the kaon is similar to 
that in the pion, and measurements of the kaon form factor, Fk{Q'^), have also recently been 
reported js^. Future measurements of Fk{Q^) at larger (~ 2 GeV^) [i^ would allow 
the first quantitative test of local Bloom-Gilman duality in strange hadrons. 
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